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TRACKING  TARGETS  WITH  SPECIFIED  SPECTRA 
USING  THE  H-PMHT  ALGORITHM 


1.  INTRODUCTION 


The  function  of  sensor-level,  multi-target  tracking  at  the  output  of  a  spatial-spectral  sensor  sig¬ 
nal  processor  (e.g.,  multi-band  remote  imaging  systems,  array  beamformers)  is  to  follow  sources, 
or  targets,  as  they  traverse  the  sensor  measurement  space.  If  target  and  noise  power  spectra  are 
specified  a  priori,  as  is  assumed  herein,  the  multi-target  tracking  problem  is  reduced  to  following 
targets  as  they  move  across  the  spatial  component  of  the  sensor  measurement  space.  In  this  case, 
target  spectra  improve  the  ability  to  estimate  target  spatial  location,  especially  in  difficult  scenarios 
such  as  crossing  targets. 

The  H-PMHT  (Histogram-Probabilistic  Multi-Hypothesis  Tracking)  algorithm  (see  Streit  [1], 
Walsh  et  al.  [2],  Streit  et  al.  [3])  is  extended  to  treat  the  case  of  targets  and  noise  having  specified 
spectra.  The  sensor  is  assumed  to  be  a  linear  system;  consequently,  target  and  noise  spectra  are 
linearly  superimposed  at  the  sensor  output.  The  histogram  model  in  the  H-PMHT  algorithm  is 
equivalent  to  linear  superposition;  that  is,  the  foundation  of  the  spectral  H-PMHT  algorithm  is  a 
model  of  linear  superposition  in  the  sensor  output.  Another  assumption  of  the  spectral  H-PMHT 
algorithm  is  that  successive  “snapshots”  of  spatial-spectral  sensor  data  are  statistically  independent 
when  conditioned  on  the  collection  of  target  states.  Conditional  independence  models  are  tradi¬ 
tional  for  sequences  of  sensor  data  obtained  from  a  measurement  process.  The  discussion  here 
assumes  that  the  spatial-spectral  distribution  of  each  target  separates  into  a  product  of  spatial  and 
spectral  factors;  however,  the  more  general  case  in  which  target  spatial  extent  is  a  deterministic 
function  of  frequency  is  treated  in  the  appendix. 

The  presentation  of  the  spectral  H-PMHT  algorithm  assumes  familiarity  with  the  technical 
details  and  notation  of  Streit  [1]  and  with  the  method  of  expectation-maximization  (EM)  as  it 
is  applied  to  multi-target  tracking  (see  Streit  and  Luginbuhl  [4]).  Closely  related  work  involving 
parameterized  spectral  mixture  models  in  a  tracking  application  includes  Luginbuhl  [5]  and  Lugin¬ 
buhl  and  Willett  [6, 7].  To  facilitate  comparisons  with  their  work,  a  detailed  algorithm  is  presented 
for  Gaussian  spectra. 
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2.  MODIFIED  SENSOR  CELL  STRUCTURE 


The  notation  of  Streit  [1]  is  adopted  here,  but  it  is  specialized  to  treat  sensor  cells  that  are 
multi-dimensional  rectangles,  i.e.,  Cartesian  products  of  spatial  and  spectral  cells.  The  sensor  cells 
C  =  {Ci, . . .  ,Cs}  are  the  Cartesian  products  of  the  U  disjoint  spatial  cells  {V\, . . . .  T>u)  and  the  V 
disjoint  spectral  cells  {£i, . . . ,  £v}-  The  particular  choice  of  spatial  and  spectral  cells  is  application 
dependent.  The  total  number  of  sensor  cells  5  =  UV,  and  every  cell  Ct  can  be  written  in  the  form 

Ci  =  Vi  x  £j  (1) 

for  some  (unique)  choice  of  the  cells  T>i  and  £j.  Let  V  =  V\  U  •  •  ■  U  I >u  =  RdlTnV  and  £  = 
£\  U  •  •  •  U  £v  =  Rdim£,  where  R  denotes  the  real  number  line. 

The  set  of  sensor  cells  from  which  measurements  are  available  at  time  t  is  denoted  by  M(t). 
This  set  may  vary  from  scan  to  scan  and  may  be  an  arbitrary  subset  of  C.  Discussing  general 
measurement  sets  M(t )  requires  cumbersome  notation;  instead,  a  special  case  is  presented,  and 
from  this  discussion  the  general  M(t)  case  is  easily  treated.  It  is  assumed  that  the  number  of  cells 
L(t)  in  M(t)  is  given  by 

L{t)  =  U{t)V(t),  1  <  U{t)  <U,  1  <  U(t)  <  V, 

and  that  M (t)  is  the  Cartesian  product  of  the  spatial  cells  {Di(f), . . . ,  Dul[t) (t)}  and  the  spectral 
cells  {Ei(t), . . . ,  .EV(t)  (£) } ,  where  Roman  fonts  are  used  (in  place  of  corresponding  script  fonts) 
to  denote  that  these  lists  are  linearly  ordered.  The  (i,  j)th  cell  in  M(t)  is  thus  Di{t)  x  Ej(t).  The 
sensor  measurement  vector  at  time  t  is  denoted  by 

%t  =  \2t,i,i>  •  •  •  1 2t,tr(t),v(t) }  > 

where  is  the  output  spectral  power  of  the  sensor  at  time  t  in  the  cell  Di[t)  x  Ej  ( t ) .  To  facilitate 
the  discussion  of  truncated  cells,  the  linear  ordering  of  the  spatial  cells  is  extended  (arbitrarily)  so 
that 


{Di(t), . . . ,  D^(t)(t),  Du(t)+ 1(*)>  •  •  •  i  Du{t)} 

is  an  ordered  list  of  all  U  cells,  and  similarly  for  the  spectral  cell  ordering.  These  orderings  change 
from  scan  to  scan  to  accommodate  different  subsets  M(t). 

The  spatial-spectral  sensor  cell  structure  in  equation  (1)  facilitates  simplifications  to  the  basic 
equations  of  Streit  [1].  Let  Xt  denote  the  set  of  target  states  at  time  t.  The  cell  probability  for  the 
(■ i,j)th  cell  takes  the  form  (cf.  [1,  equation  (6)]) 

Pij{Xt)  —  [  f{u,v\Xt)dudv,  (2) 

J  Di(t)xEj{t) 
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where  the  sample  PDF  f(u,v\Xt)  is  defined  over  all  (u,  v)  e  RdimV  x  i?dim£  =  RdimC  by  the 
mixture  (c/.  [1,  equation  (30)]) 


f(u,v\Xt)  =  ^tkGk(u,v\Xt), 

k= 0 

and  where  the  component  Gk(u,  u|Xt)  corresponds  to  target  A:  if  k  >  1  and  to  noise  if  k  =  0. 
The  total  sensor  probability  (c/.  [1,  equation  (7)])  at  time  t  becomes 

m  v(t) 

P(X<)  =  E  E 

t=l  i= 1 

The  expected  measurement  ztij  takes  the  form  (cf.  [1,  equation  (55)]) 


_  _  j  ztij ,  for  1  <  i  <  U(t),  l<j<V(t) 

Ztii  ~  \  W\  for  +  1  ^  ^  V,  V®  +  1 


<  j  <  V, 


where 


m  m 

izu-EE** 

t=i  j=i 

The  auxiliary  functions  become  (cf.  [1,  equation  (60)]) 
m  r  u  v 


v  y  __  /» 

Qt*  =  22  EE  FTri  /  Gk(u,v\x'tk)dudv  fl^logTr tk 

t=0  L“I  i=1  J Di(t)xEj(t) 


and  (cf.  [1,  equation  (61)]) 


Qkx  =  2  1OgPHt,fc|Et_1,fc(xtJb!xt_i,fc) 


T  V  V 


+  EEE  ptk(ZY'\  [  Gk(u,  v|a4)  log  Gk(u,  v\xtk)  du  dv.  (7) 

trn  1  i= 1  jZ 1  JDtWxBjW 

If  the  noise  component  model  is  specified  a  priori,  the  auxiliary  function  QkX  is  needed  only  for 
k  =  1, . . . ,  M. 
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3.  MODIFIED  AUXILIARY  FUNCTION 


The  spectral  PDF  of  target  k  is  denoted  by  Sk(v),  so  that 

J  Sk(v)  dv  =  1. 

The  spectral  PDF  is  equal  to  the  traditional  power  spectrum  normalized  so  that  its  integral  over  £  is 
one.  If  the  desired  noise  spectral  PDF  5o(u)  is  white,  it  is  constrained  to  integrate  to  one  by  making 
it  constant  over  an  appropriately  specified  finite  subset  of  £.  To  simplify  the  cuirent  discussion, 
target  spatial  and  spectral  characteristics  are  assumed  to  separate,  i.e.,  the  component  GJu 
of  the  sample  PDF  factors  into  the  form 


Gk{u,v\xtk)=Sk{v)gk(u\xtk),  k  =  l,...,M,  (8) 

where  Qk{u\xtk )  is  the  spatial  PDF  of  target  k.  The  more  general  case  in  which  the  spatial  PDF 
depends  on  frequency  is  discussed  in  the  appendix.  The  product  form  (8)  enables  multiple  integrals 
over  Di(t)  x  Ej(t)  to  be  rewritten  as  products  of  integrals,  so  that 

JdMxe  (t)  °k{'U'  V'X’tk)  dudv=  L.,  Sk^  dv  f  Gk(u\x'tk)  du  (9) 

J  Ej (£)  J £. (t)  J D. (Q 

and,  using  the  mixture  (3)  and  the  definition  (2), 


~  ^  Kk  /  Sk(v)dv  f  Qk{u\x'tk)  du.  (10) 

k—0  J  Di(t) 


Similarly, 

/  Gk{u,v\x'tk)\ogGk(u,  v\xtk)dudv 

=  /  Sk(v)dv  /  &k(‘U'\x'tk)  log  Gk(u\xtk)  du 

+  /  gk{u\x'tk)du  I  Sk(v)  log Sk{v)dv.  (11) 
JEj(t ) 

Substituting  (9)  into  (6)  gives 


M 


q** — 


U  /  V  = 


EE 


k=o  Li=i  \y=i 


Zti3  lEj(t)  $k(v)dv 


'Dili) 


Gk{u\x'tk)  du 


Kk  log  *tk,  (12) 
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and  substituting  (11)  into  (7)  gives 


QkX  =  ^2  pffij  l°g  P3t,t (xtk  ) 

,  V-  _/  V'  ('£'  2tij  k  Sk^  dv 
+  71  tk  1  z^  p-  (x') 

t=  1  i=  1  \j=l 


&(u|a4)log^(tt|xtfc)<2u.  (13) 


The  second  term  in  (1 1)  contributes  an  additional  term  to  (13),  but  it  is  omitted  because  it  depends 
on  x'tk  and  not  on  xtk  and  is  not  needed  in  the  M-step  of  the  EM  method. 


4.  LINEAR  GAUSSIAN  CASE  WITH  SPECIFIED  TARGET  SPECTRA 

■  ST?f  ^tk\  3X6  ?btained  usin§  a  LaSrangian  multiplier  technique  in  the  same  manner  as 

in  Streit  [1],  The  result  is,  from  (12),  for  k  =  0, 1, . . . ,  M, 


u  /  v 


=!ee 


^tij  )  Sk(v)  dv 

PM) 


where 


*  =  E  E 

*=0  Li=l  \j=l 
U  V 

= 

i=l  j= 1 


Zti3  §Ej{t)  Sk(y)  dv 

PM) 


Qh(u\x'tk)  du , 


[  Ok{u\x'tk), 


an  identity  obtained  by  making  the  sum  over  k  innermost  and  using  (10). 

Estimating  the  state  variables  requires  assuming  specific  parametric  forms.  Linear  Gaussian 
target  and  measurement  process  models  are  adopted  here.  These  forms  are,  for  k  =  1, ... ,  M, 

Pzt,k\2t-i,k  (xtfc|zt-i,*)  =  N(xtk\ Ft- 1,*  £t-i,fc,  Qt-iik:),  (16) 


Gk{u\xtk)  —  N{u\  Htkxtk,Rtk)-  (17) 

The  noise  component  can  be  very  general;  here,  the  noise  is  assumed  to  be  specified  a  priori  at 

^  *„  'V'1'  dr“  ”°ise  PDF  «>W*«)  3  «>(**)  is  spatially  white,  i,  is  constrained  to 
ntegrate  to  one  by  making  it  constant  over  an  appropriately  specified  finite  subset  of  V. 

Estimates  for  the  state  variables  can  be  obtained  by  setting  the  gradient  of  the  auxiliary  function 

ZZ°  and  S°Img;  however>  an  alternative  approach  is  taken  here  because  it  exploits  the 
Kalman  filter  as  an  efficient  computational  algorithm.  Completing  the  square  on  the  state  variables 

in  *  “  **  exPression  for  Qkx  and  exponentiating  the  result,  as 

m  Streit  and  Lugmbuhl  [4],  gives  the  expression  (c/.  [1,  equation  (75)]) 

T 

e  kX  K I]  Ft.!,* i(_w,  Qt_lit)  jV (4;  a;tt)  (ig) 

where  (c/.  [1,  equations  (70)  and  (71),  respectively]) 

X  P{X[+l)  ^ 

®tk  HZt+ill  ®tk’ 


7 


5  Rtk 

Rtk  =  -r—,  l  <t<T, 

77  tk  vtk 


and  the  synthetic  spatial  measurements  are  (cf.  [1,  equation  (74)]) 


u  /  v  = 


Vtk 

i=  1  \j=l 


*tij  fE.(t)  Sk(v)  dv 

PiiiX't) 


u  J\f(u\  Htk  xtk,  Rtk )  du, 


where  (cf.  [1,  equation  (72)]) 


Vtk  ~  it,  (it  ~  - Ej(t)  —— — 

i=l  \j=i  Pij(X't) 


ff(w,  Htk  xtk,  Rtk)  du. 


“  iltk}  afe’  therefore>  efficiently  computed  via  a  Kalman  filter  because  of 
he  form  of  the  likelihood  function  (18).  Further  details  are  given  in  Streit  [1]  for  converting  these 
expressions  for  a  batch  algorithm  into  a  recursive  algorithm  for  a  sliding  batch. 

(19)  “  "°W  reWritten  as  a  convex  combination  of syn- 

Stki(j)  =  ztij  —  W  dv  fp,(t)  N{u\  Rtk  x’tk,  Rtk)  du 

3  E^o  <  4(t)  SK{v)dv  JDi(t)  ff  (u\  HtK  XtK,  RtK)  du  ’  j  =  1’"->K  • 

The  broadband  power  of  the  synthetic  spectrum  Stki(j)  is  therefore 

v 

1 1 ^tki  1 1  =  ^  $tki(j),  i  =  1, . . . ,  U.  (22) 

j= 1 

Refemng  to  as  target  broadband  power  in  spatial  cell  DAt)  is  reasonable  because  it  is 

non-negative  for  all  indices  t,  k,  and  i,  and  because  S 

M  v 

X*  11^*11  = 

k= 0  j— i 

is  the  total  (expected)  power  in  D,(t).  The  target  synthetic  spectra  (21)  are  insightful  functions  that 
££££ t*„“n  *"  °f  tad»  — •  ™ey  are  also  potentially 


The  spatial  cell  centroid  ztki  for  target  k  is  given  by 


~  _  4(t)  u^(u;  Htk  x'tk,  Rtk)  du 
lDi{t)  N(u;  Htk  x'tk,  Rtk )  du 
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Each  centroid  ztki  lies  in  the  spatial  cell  Di(t)  because  it  is  the  mean  of  a  PDF  whose  support  is 
Di(t).  Multiplying  and  dividing  in  (19)  and  (20)  by  the  integral  over  Di(t),  substituting  (20)  into 
(19),  and  rewriting  the  result  gives 


ztk 


EiLi  ll^ull 


(24) 


Thus,  from  (24),  the  spatial  synthetic  measurement  ztk  for  each  target  A;  is  a  convex  combination  of 
the  spatial  cell  centroids  {ztki  :  i  =  1, and  the  coefficients  of  the  convex  combination  are 
proportional  to  the  broadband  powers  j  ||<Stfci  ||  :  i  =  1, . . . ,  £/  j  of  the  kth  target’s  synthetic  spectra. 

Writing  the  updated  mixing  proportions  (14)  in  terms  of  the  target  synthetic  spectra  gives 


or 


_  EL  g.i  SmU) 
**  EL  EL  EL 


(25) 


*tk  = 


E£=i  T.U  5*0) 
EU 


(26) 


results  that  are  easily  verified  by  direct  substitution.  Expression  (25)  is  insightful  but  computation¬ 
ally  inefficient;  however,  expression  (26)  is  potentially  useful  in  applications. 
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5.  LINEAR  GAUSSIAN  CASE  WITH  PARAMETERIZED  TARGET  SPECTRA 


In  this  section  it  is  assumed  that  a  parametric  form  of  the  target  spectral  PDF  is  specified, 
and  that  the  spectral  parameters  are  estimated  from  the  measured  data  using  the  EM  method.  Let 
Sk(v,  ytk)  denote  the  spectral  PDF  of  component  k  at  time  t,  where  ytk  is  the  spectral  parameter 
vector.  Thus,  (8)  becomes 


Gk{u,v\xtk-,ytk)  =  Gk(u\xtk)  Sk(v,ytk), 


and  (11)  becomes 


/  Gk(u,v \x'tk;  y'tk)  log  Gk{u,  v\xtk\  ytk)  du  dv 

J  Di(t)xEj(t) 

=  /  Skfay'a)  dv  /  Qk{u\x'tk)\ogQk{u\xtk)  du 
J  Ej  (t)  J  D,(t) 

+  [  Qk{u\x'tk)du  [  Sk{v,y'tk)  log Sk{v;ytk)dv.  (28) 
JDiit)  J  Ej  (t) 

Letting  Yt'  =  {y1^,  y'tl, . . . ,  y'm}  denote  the  vector  of  the  current  parameter  values  for  noise  (k  —  0) 
and  targets  ( k  >  0),  the  integral  (9)  becomes 


'DiWxEjd) 


Gk[u,v\x'tk,y'tk)dudv  =  [  Sk{v;y'tk)dv  [  gk{u\x'tk)  du,  (29) 

JEj{t)  JD^t) 


and  the  total  probability  (10)  becomes 


r  /» 

pij(XvYt)  =  ^2  ***  /  Sk{v,y'tk)dv  /  gk{u\x'tk)  du.  (30) 

fc= 0  JDi{t ) 

Using  expressions  (27)-(30)  gives  rise  to  slightly  modified  forms  of  the  H-PMHT  auxiliary  func¬ 
tions  (6)-(7).  The  auxiliary  function  for  the  mixing  proportions  is 


a.  - E  [e  (e %  Vu- y-)!t) ~ )  [  «*«•*.  (3d 

k=0  [i=l  Vi=i  *t )  )  JDi{t)  J 

and  it  leads  to  essentially  the  same  estimator  (14).  The  auxiliary  function  for  the  state  variables  is 

Qkx  =  p=-*,k  {xtk  |xt_lifc) 

JL,  Ztij  Sk(v;yL)dv\  r 

+  H  ntk  2  - P  (X>  Y'\ - I  /  GAu\x'tk)  \ogGk(u\xtk)  du  (32) 

t=l  i=l  \i=l  Yt)  J  JDiit) 


k= 0  Li=l  \j= 1 
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so  that,  for  linear  Gaussian  models  (16)— (17),  state  estimates  are  computed  via  a  Kalman  smooth¬ 
ing  filter  in  the  same  manner  as  in  Streit  [1],  The  auxiliary  function  for  the  spectral  parameters 


Qkv =^2Kk 


‘=1  3= 1  \i= 1 


*u3  /r>i(t)  Gk(u\x'tk)  du 

PiAX'tX) 


Sk(v,  y'tk)  log Sk{v,ytk)dv.  (33) 


This  expression  is  obtained  from  the  second  term  in  (28),  but  it  is  omitted  if  the  spectral  parameters 
are  not  estimated. 

The  spectral  PDF  of  target  k  is  now  assumed  to  be  Gaussian  distributed;  that  is, 

Vtk)  =N(v,  ntk,  Stk),  (34) 

where  ytk  =  {iltk)  Stk}.  The  parametric  form  of  the  spatial  dependence  function  Gk{u\xtk)  is  left 
unspecified  here.  The  synthetic  target  spectra  in  the  current  step  of  the  algorithm  are  given  by 

Stki(j)  =  Zf  ^ tk  ^ Ej ®  ^tk’  ^ik)  dv  Ip^t)  @k(u\x'tk)  du 

S-o  ^tK  S Ej (t)  N(v,  ntK,  StK)  dv  Gk (u\x'tK)  du 

If  (34)  holds,  Q^y  is  easjly  maximized  by  setting  its  gradients  with  respect  to  fj,tk  and  Stk  to  zero 

and  solving  for  $tk  and  Stk.  The  resulting  estimators  separate  hierarchically  in  the  usual  wav  and 
they  are  given  by 


^tk~x^J2Yl  StkxU) 

*  i=l  j=l 


■x  /bj(«)  vJsr{v\ y!tk,S'tk)dv 
$Ej{t)  Vtk’  $tk )  dv 


Stk  =  —  S  (i)  ^  ~  ^tk^v  ~  vtk)* M{v,  Mtfc,  S'tk)  dv 

Xtk  it  it  m  JEj{t)^’V'tk,s[k)dV 


\ tk  ~~f  .. 

1=1  j=l 


where  the  normalization  constant  is  given  by 


u  v 


Xtk=YlJ2  StkiU)- 


»=i  j= i 

These  estimators  are  readily  interpreted  as  convex  combinations  of  cell-level  mean  and  variance 
contributions. 

Explicit  spectral  parameter  estimates  can  also  be  derived  for  other  families  of  parameterizations 
of  the  target  spectral  PDFs.  For  example,  the  spectral  PDF  can  be  a  Gaussian  mixture.  For  further 

chapter^]  ^  miXtUr6S  t0  estimate  spectral  PDFs  from  histogram  data,  see  Luginbuhl  [5, 


6.  SUMMARY 


The  incorporation  of  specified  target  and  noise  power  spectra  directly  into  the  multi-target 
tracking  problem  is  accomplished  using  the  H-PMHT  algorithm.  The  synthetic  H-PMHT  mea¬ 
surements  for  this  application  are  a  pair  of  closely  related  quantities,  namely,  synthetic  spatial 
measurements  for  each  target  and  synthetic  spectra  for  each  target  in  every  spatial  cell.  These  syn¬ 
thetic  data  are  used  to  write  the  spectral  H-PMHT  algorithm  as  a  bank  of  recursively  reweighted 
smoothing  filters.  The  resulting  spectral  H-PMHT  algorithm  constitutes  a  potentially  important 
extension  of  existing  tracking  methodologies  that  may  be  useful  in  several  application  areas — for 
example,  hyperspectral  data  obtained  by  remote  imaging  systems. 
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APPENDIX 

FREQUENCY  DEPENDENT  SPATIAL  DENSITY  FUNCTIONS 

The  general  case  in  which  the  target  spatial  density  varies  with  frequency  is  discussed  in  this 
appendix.  It  is  assumed  that  the  spectral  cells  [Ei(t), ... ,  Ey(t) (i) }  composing  the  measurement 
set  A4  ( t )  are  bounded  for  all  t ,  a  physically  realistic  sensor  model.  The  factorization  (8)  is  replaced 
(using  Bayes  Theorem)  by  the  conditional  product 

Gk(u,  =  <Sfc(u)  Qk(u\v,  xtk).  (38) 

Integrating  (38)  over  the  cell  Dt(t)  x  Ej(t)  gives  (cf.  (9)) 

/  Gk(u,v\x'tk)  dudv  =  /  <<Sjfe(t;)  f  Gk(u\v,  x'tk)  du\  dv.  (39) 

The  term  in  braces  in  (39)  is  a  function  of  v  alone.  It  is  assumed  that  this  term  is  continuous  on 
Ej{t),  the  closure  of  set  Ej(t).  (This  assumption  is  reasonable  in  many  applications,  but  it  does 
have  potentially  important  implications.  For  example,  if  the  spatial  function  Gk(u\v,  x'tk)  is  a  con¬ 
tinuous  function  of  v,  the  spectrum  Sk(y)  is  continuous  and  so  must  correspond  to  a  signal  having 
no  idealized  narrowband  components.)  Applying  the  Mean  Value  Theorem  of  the  calculus  to  the 
integral  of  this  term  gives,  for  some  point  Vj  ( t )  €  Ej  (f ) , 

I  $k(v)  J  Qk{u\v,x'tk)  dudv  =\Ej(t)\  Sk(vj(t))  f  Gk{u\vj(t),x'tk)  du,  (40) 

JEiW  JDi{t)  JDi(t) 

where  |£)(f)|  =  fs^t)  <  ^  ^  spectrum  is  constant  over  Ej(t),  then 

$k(vj(t))  =  f  Sk{v)dv.  (41) 

If  the  spectrum  is  not  constant  over  Ej(t),  equation  (41)  is  only  an  approximation,  but  one  whose 

quality  improves  with  decreasing  size  of  the  cell  Ej  (t) . 

Substituting  (41)  into  (40),  and  the  result  into  (39)  gives 

Jn  t*  p  Gki~U' dudv=  [  dv  [  Qkiulvjit),  x'tk)  du.  (42) 

J  Di  (t)xEj(t)  J  Ej{t)  J  Di(t) 

The  exact  location  of  the  point  Vj(t)  e  Ej(t)  is  unknown.  Replacing  it  with  an  a  priori  specified 
point  in  Ej  ( t ) ,  denoted  by  Vj  (f ) ,  gives  the  approximation 

Id  (t)xE  (t)  dudv~  L,.  Sk ^  dv  f  Sk(u\vj{t),x'tk)  du.  (43) 

J  Di{t)xEj(t)  J Ej{t)  J Di(t) 
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u  kno,wn’ the  aPProximation  (43)  may  be  used  to  replace  the  factorization  of  equa¬ 
tion  (9)  throughout  the  discussion.  4 

,s„The,?Irary  ‘Tpact  °f  “Sing  the  approximation  (43)  on  the  spectral  H-PMHT  algorithm  is 

rim  spatral  lntegra]s  now  depend  on  ^  specmd  tadex  .  vja  ^  parameKr  -  {()_  so 

double  summations  cannot  be  evaluated  as  efficiently  as  in  the  frequency-independent  case  The 
numerical  details  depend  on  the  parameterization*.  For  example,  if  the  specified  spa  aT form  of 
yk{u\vj{t),xtk)  is  given  by  (cf  equation  (17))  P 

Gk ( u\vj (t) ,  xtk)  =  M  (u-  Htk  xtk ,  Rtk/  (7  vj (t) )2) ,  (44) 

where  %  (t )  is  the  midpoint  of  the  cell  Ei(t)  and  7  is  a  specified  constant,  then  foe  estimate  of  the 
synthetic  target  spectrum  Stki(j)  becomes  (cf  equation  (21)),  for  j  =  1, . . . ,  V, 

Stki(j )  =  ztij  ~  ImL  Sk^dv  Spdt)N  (u\Htkx'tk,  (yvj(t))~2 Rtk)  du 

T,k=o  Kk  fEj(t)  s*(v)  dv  fDi(t)  Af  (U;  Hu  x[K,  (jVj(t))~2  Rtk)  du  (45) 

detoilsTeomklIdSUrementS  %  ^  m°dified  by  the  P^etenzation  (44);  however,  further 
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